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Abstract 

We calculate the top-quark-induced three-loop corrections of 0{agGpMf) to the 
Yukawa couplings of the first five quark flavours in the framework of the minimal 
standard model with an intermediate-mass Higgs boson, with mass Mjj <C 2Mt- 
The calculation is performed using an effective-Lagrangian approach implemented 
with the hard-mass procedure. As an application, we derive the 0(o^G_pM t 2 ) correc- 
tions to the H — > qq partial decay widths, including the case q = b. The couplings 
of the Higgs boson to pairs of leptons and intermediate bosons being known to 
0{a 2 s GpMf), this completes the knowledge of such corrections in the Higgs sector. 
We express the results both in the MS and on-shell schemes of mass renormaliza- 
tion. We recover the notion that the QCD perturbation expansions exhibit a worse 
convergence behaviour in the on-shell scheme than they do in the MS scheme. 
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1 Introduction 



One of the longstanding questions of elementary particle physics is whether nature makes 
use of the Higgs mechanism of spontaneous symmetry breaking to endow the particles 
with their masses. The Higgs boson, H, is the missing link sought to verify this the- 
oretical conjecture in the standard model (SM). The possible range of the Higgs-boson 
mass, Mff, is constrained from below both experimentally and theoretically. The failure 
of experiments at the CERN Large Electron-Positron Collider (LEP 1) to observe the 
decay Z — > ffH has ruled out the mass range Mh < 65.6 GeV at the 95% confidence 
level 0. Depending on the precise value of the top-quark mass, M t , the requirement that 
the vacuum be the true ground state provides an even more stringent theoretical lower 
bound 0. However, this bound may be somewhat relaxed by taking into account the 
possibility that the physical minimum of the effective potential might be metastable || . 
Other theoretical arguments bound Mh from above. The requirements that partial-wave 
unitarity in intermediate-boson scattering at high energies be satisfied @ or that pertur- 
bation theory in the Higgs sector be meaningful || establish an upper bound on Mh at 

about rs 1 TeV, where Gp denotes Fermi's constant, in a weakly inter- 

acting SM. The triviality bounds, i.e., Mh upper bounds derived through perturbative 
H or lattice computations by requiring that the running Higgs self-coupling, A(/x), 
stay finite for renormalization scales fi < A, where A is the cutoff beyond which new 
physics operates, are somewhat stronger. In the following, we shall focus our attention 
on the lower end of the allowed Mh range, where Mh <C 2M t , which will be accessed by 
colliding-beam experiments in the near future. 

A Higgs boson with Mh ^ 135 GeV decays dominantly to bb pairs ||. This decay 
mode will be of prime importance for Higgs-boson searches at LEP 2 ||, the Fermilab 
Tevatron [IIJ or a possible 4- TeV upgrade thereof ]TT], a next-generation e + e~ linear 
collider |f[2]| , and a future collider [13|. Techniques for the measurement of the 

H — > bb branching fraction at a y/s = 500 GeV e + e~ linear collider have been elaborated 



in Ref. [fL4|| . The branching ratio of H — > cc is roughly 30 times smaller than that of 

Once a novel scalar particle is discovered, it will be crucial to decide if it is the very 
Higgs boson of the SM or if it lives in some more extended Higgs sector. To that end, 
precise knowledge of the SM predictions will be mandatory, i.e., quantum corrections must 
be taken into account. The present knowledge of quantum corrections to the H — > qq 
partial decay widths has recently been reviewed in Ref. ||15|| . At one loop, the electroweak 
fi~6| , [P7H and quantum-chromodynamical (QCD) (T^| corrections are known for arbitrary 
masses. The leading high-M^ term, of 0{GpM H ), was first derived by Veltman |Tj| . 
In the limit Mh <C 2M t , which we are interested in here, the terms of 0(X t ), where 
X t = (G F M?/8n 2 V2),tend to be dominant. They arise in part from the renormalizations 
of the Higgs wave function and vacuum expectation value, which are independent of the 
quark flavour q [2D|. In the case of q = b, there is an additional non- universal 0(X t ) 



contribution ITq, 17, 21], which partly cancels the flavour-independent one. At two loops, 



the universal \22\ and bottom-specific p3|, |24| 0(a s X t ) terms are available. Furthermore, 
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the first |25j and second []26| terms of the expansion in M 2 / of the 0(a 2 ) five-flavour 
QCD correction have been found. As for the top-quark-induced correction in 0(a 2 ), the 
full M t dependence of the non-singlet (double-bubble) contribution |27| as well as the 
first four terms of the Mfj/M 2 expansion of the singlet (double-triangle) contribution 
have been computed. At three loops, the 0(a^) non-singlet correction is known in the 
massless approximation [f29f| . 

In this paper, we shall take the next step, to three loops including virtual top-quark 
effects. Specifically, we shall evaluate the 0(a 2 X t ) corrections to the decay widths of 
H — > qq, and in particular of H — > bb. These corrections may be divided into three 
classes, which are separately finite and gauge independent: (i) The universal correction 
originates in the renormalizations of the Higgs wave function and vacuum expectation 
value and occurs as a building block in the computation of quantum corrections to any 
Higgs-boson production or decay process. It is related to the Higgs- and W^-boson self- 
energies. In the case of the leptonic decay, H — > l + l~ , this is the only source of 0(a 2 s X t ) 
corrections |p0[ . (ii) The quark- specific correction arises from Feynman diagrams where 



the top quark only appears in a closed loop which is connected with the external quark line 
by two gluons (and, in some cases, one additional weak neutral boson). It just depends 
on the third component of weak isospin of the considered quark. (Hi) The bottom-specific 
correction emerges from the one-loop seed diagram where one charged boson is emitted 
and re-absorbed from the external bottom-quark line, by appropriately adding gluon lines 
(and, in some cases, one quark or ghost loop). It is the natural extension of the bottom- 
specific 0(a s X t ) correction to the H — > bb decay width p3|, |4j] by one order of a s . 

The next-to-leading-order QCD corrections to the top-quark-induced shifts in the 
l + l~H, W + W~H, and ZZH couplings have been found in Ref. fl3"0f . The present paper 
completes the knowledge of the 0(a 2 X t ) corrections to the SM Yukawa couplings, still 
excluding that of the top quark. The Yukawa couplings of the first four quark flavours 
receive corrections from the same class of diagrams, while, in the case of the bottom 
Yukawa coupling, an additional class of diagrams must be included. 

In Ref. [0, it was noticed that the QCD perturbation expansions of the l + l~H, 
W + W~H, and ZZH couplings and the electroweak p parameter, for which the 0(a 2 X t ) 
correction is also known [|3~H| , exhibit striking similarities. If the top-quark mass is renor- 
malized according to the on-shell scheme, the coefficients of (a s /7r)X t and (a s /ir) 2 X t are 
in each case negative and of increasing magnitude. On the other hand, the corresponding 
coefficients in the MS scheme are much smaller and of variant sign. This gave support 
to the notion that the use of the pole mass deteriorates the convergence properties of the 
QCD perturbation series, which may also be motivated from the study of renormalons 
|32|| . It is interesting to find out whether this observation is substantiated by the analysis 



of the light-quark and bottom Yukawa couplings. We shall return to this issue in Sect. [g. 

Our key result for the H — > bb decay width to 0(a 2 X t ) has recently been presented 
without derivation in a brief note This paper provides the full details of our analysis 
and also deals with the H — > qq decays, where q — u,d, s, c. It is organized as follows. 
After introducing our notations in Sect. |2|, we shall derive, in Sect. [|, a heavy-top-quark 
effective Lagrangian with two coefficient functions to be determined by diagrammatic 
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calculation. From this Lagrangian, we shall obtain a generic formula for the H — > qq 
decay widths valid through 0(a 2 X t ). In Sect. f|, we shall calculate, at three loops, the 
coefficient functions relevant for the first four quark flavours. In Sect. |5|, we shall extend 
this analysis to also include the case q = b. Our conclusions will be summarized in Sect. |]. 



2 Notations 

In this section, the notation is fixed and useful formulae, which will be necessary to numer- 
ically evaluate our results, are provided. The calculation is performed in the framework of 
dimensional regularization with space-time dimension n = 4 — 2s. The QCD gauge group 
is taken to be SU(iV c ), with N c arbitrary. The colour factors corresponding to the Casimir 
operators of the fundamental and adjoint representations are C F = (N 2 — 1)/(2N C ) and 
Ca = N c , respectively. For the numerical evaluation we set N c = 3. The trace normaliza- 
tion of the fundamental representation is T = 1/2. The number of active quark flavours 
is denoted by rif. Unless otherwise stated, we work in the MS scheme, with fi being the 
renormalization scale. 

The \x dependence of the strong coupling constant, a s (/z), and any quark mass, 
is governed by the renormalization-group (RG) equations, 
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are the first few coefficients of the QCD (3 and 7 functions. The next-to-leading-order 
(two-loop) solution of Eq. ([]]) reads 
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where A is the asymptotic scale parameter. 

The relation between the MS mass m{fj) and the on-shell mass M t of the top quark 
is given by 
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where n/ = 6, X t = \GpMj* /87v 2 y/2) , and £ is Riemann's zeta function, with values 
C(2) = tt 2 /6 and C(3) « 1.202 057. The C(a 2 ) and 0(X t ) corrections were calculated in 
Refs. |Q and |35], |36|, respectively. Iterating Eq. (||), one obtains the scale-independent 
MS mass \i t = m t (fi t ) as 
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+ CaCf 



1111 



384 



+ ic(2)(l-31n2) + |c(3) 



+ C F Tn f ( — + -C(2) 



0(a 3 s ,a s X t ,X?). (7) 



In the MS scheme, the relation between the values of a s (fi) for rif = 5 and rif = 6 is given 
by 



«i 5) (^) = afV)<!l + 



(6), 



7T 



— - In At + a; t f —\ + In At 



mi 
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where x t ([i) = Gptn 2 (yu)/87r 2 v / 2 . The constant term of 0{a 2 s x t ) in Eq. (§) represents a 
new result, which follows from the analysis of the next section. Inserting the 0(X t ) term 
of Eq. (H) into Eq. (jST), one obtains 



«i 5) ^) = arv)<!i + 



(6), 



T 



7T 
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+ O(a^) . (9) 



3 Effective Lagrangian 

In this section, we construct an effective Lagrangian for the interaction between an 
intermediate-mass Higgs boson and a quark- ant iquark pair. Therefore, in addition to 
the pure QCD Lagrangian, the couplings of the quarks to the Higgs boson, H, the neutral 
Goldstone boson, x, an d the charged Goldstone bosons, ± , must be taken into account. 
This will produce corrections proportional to X t . In this paper, we are not interested in 
corrections of 0(X 2 ), and our formulae will not in general be valid in this order. 
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As a starting point, we consider the bare Yukawa Lagrangian, 



o ' 
v' 



(10) 



where v is the Higgs vacuum-expectation value, the superscript labels bare quantities, 
and the operator J is defined as 



J 



J2m° q( fq + m t ft . 



Here, q runs over u, d, s, c, and b. It is easy to see that J is a finite operator, in the 
sense that no additional renormalization constant is needed. Our aim is to construct 
the equivalent expression for J in the effective theory where the top quark is integrated 
out. Because only the leading terms in M t are considered, Ly may be written as a linear 
combination of three physical operators with mass dimension four, namely, 
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where 
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0/ 



-,0' 



are bare operators in the rif = 5 effective theory and C° 



(13) 

fi) (i = 1, 2q, 3q) are 

their bare coefficient functions. Here, G aflv is the colour field strength, is the covariant 
derivative, and the parameters and fields of the effective theory are marked by a prime. 
Notice that 0' 3q vanishes by the fermionic equation of motion and may be omitted once 
and C 2q are determined. 

The relations between the parameters and fields in the full and effective theories read 



,0/ 
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G 

m 
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G 



m° 

m,q" "q 




(14) 



where G af j, denotes the gluon field. The renormalization constants ( 2q i Cm, 9 ) an d C3 pl & Y 
an important role in the determination of the coefficient functions Cf. As may be seen 
using the method of projectors P7[ , they may be calculated from the vector and scalar 
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parts of the quark self-energy, Ey(p 2 ) and £s(p 2 ), and the transverse part of the gluon 
self-energy, H(p 2 ), via 

C 2 ° 9 = l + ^(0), 
1 - ££(0) 



m,q 



1 + 5^(0)' 

i + n°*(o), 



(15) 
(16) 
(17) 



where the superscript t indicates that only diagrams containing the top quark have to 
be considered and the superscript reminds us of that the parameters are still in their 
bare forms. In our convention, the bare quark and gluon propagators are proportional 

to [y (1 + Y? v (p 2 ) + isY? A {j) 2 )) - m° q (1 - S^p 2 ))]" 1 and [p 2 (l + U°(p 2 ))]-\ respectively 
Notice that the axial- vector part of the quark self-energy, E^p 2 ), does not enter our 
analysis because it leads to corrections of 0(x 2 ). The renormalization constant (® only 

enters the stage via Eqs. (|8]) and (^|). Through 0(a 2 x t ), we have (® = ((®) . 

By means of the higher-order formulation [|38| of a well-known low-energy theorem 
(LET) E{| in connection with the method of projectors [57|, one may derive the following 
relations which allow one to calculate the coefficient functions: 



S2qSm,q \y\ 
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dm 
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(19) 
(20) 



Here it is understood that the operator (m°d/dm®) only acts on those appearances of 
m° which remain after every l/v° is saturated by one power of m° to give m®/v°. As 



may be seen from Eqs. (|II| )-(f20|), the calculation of the coefficient functions involves the 
quark- and gluon-propagator diagrams containing the top quark, with nullified external 
momenta. By means of the LET, the corresponding vertex diagrams may be generated 
by attaching an external Higgs-boson line with zero momentum to the quark lines. To the 
order of interest here, this requires the evaluation of one-, two-, and three-loop tadpole 
diagrams. They are calculated with the help of the program package MATAD, written in 
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FORM P0[ , which makes use of integration- by-parts identities developed in Ref. 

So far, we have constructed an expression for J in terms of bare operators and coeffi- 
cient functions, which both still contain poles in e. As is well known, different operators of 
the same dimension and quantum numbers in general mix under renormalization. Specif- 
ically, the renormalized operators, which will be denoted by square brackets, are related 



to the unrenormalized ones according to |42 



[0[] 
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(21) 
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where Z g and Z m are the coupling and mass renormalization constants, respectively, in 
pure QCD with n/ = 5 active flavours. It hence follows that 



fj, 2 d 
dfi 2 

fi 2 d 
djj? 



[0[] = - 



da. a. 



a s d 



™ +4 ^ 



E K 



o' 



2q 
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Note that 0' 3 does not mix with 0[ and 0' 2q . On the other hand, the coefficient functions 
are renormalized according to 



Ci 



1 



l + 2(a s d/da s )\nZ g 

4(a s d/da s )\nZ m 
l + 2(a s d/da s )\nZ g 
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Consequently, the n/ = 5 effective Lagrangian ( |12|) takes the form 



Ci [o[] + £ c 29 + Cs, 



3</ 



(24) 



The new coefficient functions and the operators are individually finite, but, with the 



exception of 



O' 



2q 



O' 



3g 



will be 
H — > qq 



they are not separately RG invariant. From now on, 
omitted because it vanishes on mass shell and thus does not contribute to the 
decay widths. 

It is desirable to arrange for the renormalized coefficient functions and operators to 
be separately fi independent up to higher orders. This may be achieved by reshuffling 



the terms in Eq. (|24j). Exploiting the RG invariance of 
energy-momentum tensor, 



O' 



2o< 5 > [ lJ 



1 



and of the trace of the 



(25) 



where all quantities are defined in the rif = 5 effective theory, we may construct two new 
operators which are indeed RG invariant [43], e.g., 

2tt fn{3^ 



O' 



[Oil " 27^ 5) E K 



O' = O' 



'2q 



Thus, the relevant part of the nt — 5 effective Lagrangian may be written as 
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where 



a 
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(5) 
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As the new coefficient functions and operators are separately RG invariant, we may choose 



u,= u, t for C g and C q and /z = Mh for 



and 



0' 



cff 



In order to calculate the H — *■ gg decay width, it is more convenient to re-express £|J 



in terms of the operators [O^J and 



O' 



2g 



At the same time, the separation of the scales 



Mh and u. t is kept, so that Eq. fl2"T|) becomes 
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(29) 



where 



Ci(nt,M H ) 



C 2q (» t , M H ) = ( 7 ^(^) - 7^ 5) (Mh)) CxCAi*) + C 2q (» t ). 



(30) 



Notice that, for our purposes, C\ and C 2q are needed up to 0(a s x t ) and C^afa^), respec- 
tively. It is instructive to expand Eq. (|3(]) in terms of a^(fj, t ), keeping only terms relevant 
for our 0(a 2 s Xt) calculation. Observing that C% starts at 0(a s ), while C 2q = 1 to lowest 
order, this yields 



C 1 (ih,M h ) = C 1 (jh), 

o «i 6) (^) 



C 2q (fi t ,M H ) 



71 
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The coefficients C\ and C 2q must be computed diagrammatically. This will be done 
for the light-quark flavours in Sect. |3] and for q = b in Sect. ^. In the remainder of this 
section, we present a generic formula, derived from Eq. fl2~9|) , for the H — > gg decay width, 
appropriate for M# <C 2Mj, which accommodates all presently known corrections. It 
reads 



T(H^qq)=T^r (l + 5 u 
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(33) 



is the Born result including the full mass dependence. As is well known [Tj|, we may 
avoid the appearance of large logarithms of the type ln(M|/m 2 ) in the QCD correction, 
A® CD , by taking m q in Eq. (|33|) to be the MS mass for n/ = 5, m^(/i), with fx of order 
Mn. Consequently, we may put m q = in A^ CD . We may proceed similarly with the 
quantum-electrodynamical (QED) correction, A^ ED , which then takes the form 



A QED = "Mg* 



17 3, fx 2 



(34) 



where a(fx) is the MS fine-structure constant and Q q is the fractional quark charge. In 
turn, m q b \fx) is then also shifted by a QED correction []35] from the pole mass, M q . 

denotes the weak correction with the leading 0(x t ) term stripped off. If we 
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where (? w = 1 — s 2 = M^/M|, M z is the Z-boson mass, and M w is the H^-boson mass. 
A^ CD is the well-known QCD correction in the nf = 5 effective theory [25 
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This correction originates from the class of diagrams where the Higgs boson directly 
couples to the final-state qq pair, i.e., it does not comprise the double-triangle topologies 
of Ref. [^8[[ . The latter will be discussed below. The scale fx in A^ CD must be identified 
with that of m q b \fx) in r^ rn . It is natural to choose fx = Mu in order to suppress the 
logarithms of RG origin. 

The corrections of 0(a™x t ), with n = 0, 1, 2, as well as the rif = 6 QCD corrections of 
0(a s ) and 0(a 2 ) are all contained in 5 U , Ci, and C 2 g- contains the universal 0{a™x t ) 
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corrections, which originate from the renormalizations of the Higgs wave function and 
vacuum expectation value. Specifically, we have 



fL- = 2^GfH ( 1 + S u 



(37) 



with 30 
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(38) 



where 5 2 « 0.260 434, D 3 » -3.027 009, and B A « -1.762 800 are mathematical constants 
related to certain three-loop tadpole diagrams. The formula for iV c arbitrary may be found 
in Ref. ||3(J|| . The choice fi = fit eliminates the logarithms in Eq. (|3"8"D and may thus be 
considered natural. 

As is well known |28 |, starting at (9(a 2 ), Y (H —>■ qq) also receives leading contributions 
from the qq and qqg cuts of the double-triangle diagrams where the top quark circulates 
in one of the triangles. The additional exchange of a virtual Higgs or Goldstone boson 
within the top-quark triangle of the double-triangle seed diagram gives rise to a 0(a 2 s xt) 
correction, which we must include in our analysis. In the framework of the rif = 5 
effective theory, where the top quark only appears in the coefficient functions, this class 
of contributions is generated by the interference diagram of the operators [0[] and 0' 2q 
depicted in Fig. [TJ. The absorptive part of this diagram also includes a contribution from 
the gg cut, which is well known and must be subtracted in order to obtain the desired 
0(a 2 s Xt) correction to Y (H — > qq). We so obtain 
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While S^ CD is obviously not RG invariant, the physical observable Y (H — > qq) is because, 
in Eq. ([32]), the fi dependence of S^ CD is compensated by that of C2 q , as will become 
apparent in the next section. 

As a by-product, we may derive from the n/ = 5 effective Lagrangian ( |2"9"| ) a formula 
for the H —>■ gg decay width which includes the 0(x t ) correction. Our result is 



Y(H^gg) = A gg (l + 5 u f (d) 2 , 



where 
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(40) 
(41) 
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Figure 1: Interference diagram of the operators [0[] and 
circles, which generates T™^. 



O' 



2q 



represented by the solid 



Notice that, to the order of our calculation, C\ is independent of the quark flavour. 

At this point, we should mention that the QCD correction H3] to T(H — > gg) 
also includes contributions due to qqg final states. These may also be interpreted as 
0{a i s M'j { /m 2 q ) corrections to the respective H — > qq decay widths Specifically, 

these corrections would appear in form of a term proportional to (a^/n) (M^/m 2 ) (Ci) 2 
within the square brackets of Eq. fl32|). In the following, this term will not be taken into 
account, since it is formally one order of a s beyond our considerations. 

In the next section, we shall derive C\ as well as C2 q for q ^ b. This will enable us 
to calculate F (H — > qq) through 0(a 2 x t ) from Eq. (|32|) and T(H — > gg) through 0(x t ) 
from Eq. (flOp. The case q = b is more complicated and will be treated in Sect. [S[ 



4 H —> qq decay to 0(a 2 s x t 

In this section, we calculate the 0(a 2 x t ) corrections to the H — ► qq decay widths, where 
q 7^ b. The relevant types of diagrams are depicted in Figs, ^(a)-(c). The diagrams of 
types (b) and (c) emerge from the pure QCD diagram in Fig. |2](a) by allowing for the 
exchange of one virtual Higgs or Goldstone boson. Although the pure QCD diagram only 
contributes in 0{a 2 s ) 1 it is needed for the renormalization. There, the top-quark mass 



appears logarithmically and has to be replaced according to [35, pq 



m° t = m t (l + ^x t ^j . (42) 

The counterterm thus obtained cancels the ultraviolet sub divergences of the type-(b) 
diagrams; the remaining overall divergences are removed if the LET is applied. We 
checked that this is true separately for the contributions due to the H, x, and ^ bosons. 

The diagrams of type (c), where the scalar-boson lines are attached to both the top- 
quark loop and the light-quark line, are also separately finite upon application of the 
LET. Due to the appearance of one power of the light-quark Yukawa coupling, in our 
approximation, they only contribute to the scalar part of light-quark self-energy, Ss(0). 
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Figure 2: Typical diagrams generating universal 0(a 2 Xt) corrections to V (H —>■ qq). The 
dashed line in diagram (b) represents a H, x, or (p^ boson, while in diagram (c) only the 
H and x bosons appear. In total, there are 12 diagrams of type (b) and 12 diagrams of 
type (c). 



A potential problem arises in the diagrams of type (c) if the scalar particle is the neutral 
Goldstone boson x, since then only one 75 matrix appears in each quark line. However, 
the arbitrariness of the 75 definition only affects the finite terms of the original diagrams, 
which are removed by taking the derivative according to the LET. 
In the MS scheme, the coefficient functions are found to be 



ai 6 )(u) / 1 1 



71 



I 12 



1, H 2 

- m 

2 m? 



x t 



--18C(3) + /,(6 + 36C(3)) 



(43) 
(44) 



where I q is the third component of weak isospin, i.e., I q = 1/2 for up-type quark flavours 
and I q = —1/2 for down- type quark flavours. The dependence on I q stems from the qqx 
coupling, which appears linearly in the x-exchange diagrams of type (c). There are several 
checks for Eqs. fl4*3|) and ([44]). If we ignore the RG improvement of Eq. ( |3~0D and insert 
the 0(a s ) term of C\ and the 0(a 2 ) term of Ci q into Eq. (p2|), then we recover the 0(a 2 s ) 
double-triangle contribution to Y (H — > qq) found in Ref. [p8]. The 0(a s x t ) term of C\ 
may be deduced from Ref. [ 4^| . The 0(a 2 Xt) term of C2 is new. In compliance with RG 
invariance, it does not contain a logarithm of the type \n(n 2 /m 2 ). 

Since the measured top-quark mass corresponds to the pole mass M t , it is convenient 
to directly express the perturbative expansions in terms of M t . In Eqs. (^) and (pUp, 
the radiative corrections are arranged in such a way that, in the limit Mu -C 2M t , all 
dependence on the top-quark mass is concentrated in <5 U , C 1; and C 2q . Consequently, we 
just need to substitute Eq. (|j) into Eqs. ([38|), ®), and (flip. In the case of C-iq-, this leads 
to 



^os 

U 2q 



1 + 



a 



(6) 



7T 



iln^ 
2 



M 2 
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T 



18C(3) + /,(6 + 36C(3)) 
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3 a* 2 



(45) 



where X t is defined below Eq. (^), while C° s emerges from C\ by merely replacing x t 
with X t . As for 5 U , its counterpart in the on-shell scheme, 5 n , may be found in Ref. 
We now turn to the H — > gg decay width. From Eq. ( pUp we obtain 



r(# - <w) 



.4 



yy 



a 



(6) 



(/*)' 



144 



1 + X 



7T 



t; 5 



(46) 



which is in agreement with Ref. [f5|. 



5 i7 — > bb decay to 0{a 2 s x t ) 

In this section, we upgrade the 0(a 2 x t ) calculation of Y (H — > qq) to include the case 
q = b. In addition to the diagrams in Fig. ^ we must now consider diagrams where the 
primary bottom quark branches into a top quark and a charged Goldstone boson <p ± . 
Typical examples are depicted in Fig. |]. At the three-loop level, there is a total of 54 
such diagrams. These diagrams only affect C^, while C\ remains unchanged. Because 
the top-quark-induced corrections already start at 0(xt), the renormalization constants 
and Cm, g °f Eqs. ( |i~5|) and QTSj) , respectively, do contribute to Eqs. (|19|) and (pD|), from 
which we extract C25. Thereby, it is important to observe that the interference of 0(x t ) 
and 0(a 2 ) terms contribute in 0(a 2 x t )- 



^ Xl,# y %q> 

Figure 3: Typical diagrams generating non-universal (9(a™Xt) (n = 0, 1, 2) corrections to 
r (h — > bb^j. The dashed lines represent the charged Goldstone boson In total, there 
are 4 (54) irreducible two-loop (three-loop) diagrams. 



In the MS scheme, C2& is found to be 

//■_•> \ 2 

441 



C, b = Cm + < -3 + ^M Cf (4 - £ in 4) + f ^ ' 

IX \ 4 2 77^ / \ 7T 



<3 



64 
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w x 99 w x 27, /i 2 27, 2 u 2N 
+ 9C 2 + — C 3 - — In ^ - — In 2 ^ 
4 16 8 



87 M 2 33 2 /i 2 \ /241 p a , 3 , 2 A* 2 V 

8 16 m^y y 48 mi A mi J 



3253 



aW(A*) 



-7-61n 



«i 6) (^)' 
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7T 



m, 



7T 



X 



43 

y 



69 /i 2 45 , , /i 2 



+ 30C(2) + — C(3) - — ln ^ - — ln z 



48 v 7 2 ' v 7 2 m 2 4 

where n/ = 6. Replacing m 4 by M t via Eq. (|6|) leads to 



m\ 



COS /^OS I v 
2b — °2d + A * 



_ 3 + q^M ^a + 

IX 4 \ 7T 



Cf T (-^ + 9C(2)) 



(47) 



+ C 2 F 



f + 9 C (2)g-21n2 



117 



C(3) 



-| + C (2)g + 91n2 



> + i in w) + ^te- 5c(2) -i in € 



C 2 ° d s + X,<!-3 + 
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qj 6) M , M 6) (/^ 
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71 
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(48 



To simplify Eqs. (f4"T|) and fl4"8|), we set N c = 3, insert the numerical values of the math- 
ematical constants, and eliminate the RG logarithms by choosing \x = \xt and \x = M t , 
respectively. We so obtain 



a 



2b 



l 5 ( af\y t ) 

18 I 7T 



-3-7 



7Y 



+ 



7T 



x (-28.018 + 1.154 nf) 



(49) 



^os 

°2b 



5 cW) 

18 I 7T J 



-3 + 



+ 



7T 



7T 



x (64.223 - 5.094 



(50) 



where we have displayed the coefficient of n/ = 6, for reasons which will become clear in 
a moment. 

Equation (|50|) extends the non-universal correction 5 nu defined by Eq. (12) of Ref. 
to 0{a 2 s X t ). From Eq. @ we read off that the leading 0(X t) term receives the QCD 
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correction factor 



We thus recover a pat- 



- 0.333 af\M t )/ix - 11.219 (af\M t )/^) 

tern similar to the electroweak parameter Ap |3T| and the corrections 5 U , Swwh, and Szzh 
30| to the l + l~H, W + W~H, and ZZH vertices, respectively. In fact, the correspond- 
ing QCD expansions in a < f > (Mt)/n of these four observables all have negative coefficients 
which dramatically increase in magnitude as one passes from two to three loops |3(J. On 
the other hand, if the top-quark mass is renormalized in the MS scheme at scale p = p t , 
then the respective QCD expansions in af^n^/n are found to have coefficients which 
have variant signs and nicely group themselves around zero |30| . In the case of C 2 b, the 



QCD correction factor reads 



We conclude 



1 + 2.333 af\u, t )/n + 7.032 (a® 

that, also in the case of the bbH interaction, the QCD expansion in the on-shell scheme 
exhibits a worse convergence behaviour than the one in the MS scheme. However, the 
difference is less striking than in the previous four cases. 

Finally, we would like to test Broadhurst's rule concerning the naive non-abelianization 
of QCD [|47|. Guided by the observation that the ^/-independent term of Pq in Eq. (|3]) 
emerges from the coefficient of rif by multiplication with —33/2, Broadhurst conjectured 
that this very relation between the ri/-independent term and the coefficient of rif approx- 
imately holds for any observable at next-to-leading order in QCD. In Ref. [ 30fl , this rule 



was applied to Ap, 5 U , Swwh, and Szzh, and it was found that, in all four cases, the signs 
and orders of magnitude of the n^-independent terms are correctly predicted. Except for 
Szzh, these predictions come, in fact, very close to the true values. If we multiply the 
coefficients of rif in Eqs. (|49|) and (|50"D with —33/2, we obtain —19.041 and 84.055, which 
has to be compared with the respective n/-independent terms, —28.018 and 64.223. Once 
again, the signs and orders of magnitude of the n^-independent terms, which are usually 
much harder to compute than the coefficients of rif, come out correctly. 



6 Discussion and summary 

In this paper, we calculated the three-loop 0(a 2 s X t ) corrections to the H — > qq decay 
widths of the SM Higgs boson with mass Mh <C 2M t , including the case q = b. To 
this end, we constructed a rif = 5 effective Yukawa Lagrangian by integrating out the 
top quark. This Lagrangian is a linear combination of dimension-four operators acting in 
QCD with rif = 5 quark flavours, while all M t dependence is contained in the coefficient 
functions. We renormalized this Lagrangian and, by exploiting the RG invariance of the 
energy-momentum tensor, rearranged it in such a way that the renormalized operators 
and coefficient functions are separately independent of the renormalization scale, p, to 
the order considered. The RG-improved formulation thus obtained provides a natural 
separation of the rif = 5 QCD corrections at scale p = M H and the top-quark-induced 
rif = 6 corrections at scale p = M t , in the sense that the final result does not contain 
logarithms of the type ln(M^/M H ) if the rif = 5 and rif = 6 corrections are expanded in 
af ) (M H ) and af\M t ), respectively. 



In contrast to the two-loop 0(a s X t ) case |23[| , where it was sufficient to consider just 
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one term in the Lagrangian, we needed to take into account three types of operators 
and to allow for them to mix under renormalization. The mixing terms are related to 
the 0(a 2 ) double-triangle contribution considered in Ref. [28] and extend the latter to 
0(a 2 s X t ). 



Similarly to Ref. |23[ , we could take advantage of a low-energy theorem to simplify the 
calculation of the coefficient functions. This allowed us to relate a huge number of three- 
loop three-point diagrams to a manageable number of three-loop two-point diagrams. 
Specifically, we had to compute 24 irreducible three-loop two-point diagrams for q ^ b 
and, in addition, 54 ones for q = b. Such a theorem is not available for the gauge 
interactions, which might explain why three-loop 0(a 2 X t ) corrections have not yet been 
calculated for the Z — > qq decay widths, including the important case of Z — > bb, which 
has recently attracted much attention in connection with the so-called anomaly. 

To illustrate the effect of the 0(a 2 X t ) corrections to the H — > qq decay widths, we 
rewrite Eq. (B2I) as 



r (h -> qq) = r*™ 



where 



(l + Af D ) (1 + A- eak [ t=o J (l + AQ CD ) (l + Aj) + ~ 9 QCD H 



(51) 



A* = (l + 5 U ) (C 2q ) 2 -1 



l + U C X C 



-2q 



(52) 



contain the leading M t dependence and refer to the n f = 6 theory, while A^ CD and H^ CD 
are confined to pure QCD with n/ = 5. We now collect our final results for A* and 
S*. For simplicity, we undo the RG improvement of Eq. (|3~0"D and employ the expanded 
versions of C\ and C 2q given in Eq. ([H]) instead. In the MS scheme with \i — fi t , we then 
have 
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(53) 



where q = u,d, s, c. Clearly, A* = A^ and A* = A^. Introducing the top-quark pole mass 
M t and choosing \i — M t , we find the equivalent expressions 
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t,os 
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A i,OS 



/ 8 M, 2 N 
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' af\M t ) \ 2 (5 2 

U + 3 ln M|. 
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7T 
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( 2, Mr 
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(54) 



The expressions for S*' os and S^' os emerge from those of S* and H£ by merely replacing 
fit and x t (fj,t) with M t and X 4 , respectively. 

At this point, we wish to emphasize that our calculation is only valid in the limit 
Mh <C 2M t , where the LET is applicable. We expect that this is a good approximation 
for a Higgs boson in the intermediate-mass range; but, strictly speaking, this is just a 
conjecture, which can only be proven once the calculation of the full mass dependence 
is available. However, we can check the validity of this assumption at one loop, where 
the radiative corrections are fully known [|Tj], |T7). Assuming M t = 175 GeV and Mh = 

amounts to 81% of 



weak 



x t =0 



100 GeV, we find that the 0(X t ) term subtracted from A 

the top-quark- induced correction inherent in A™ eak . It is reasonable to expect that this 
feature carries over to higher orders. 

Notice that the branching ratios of the Higgs boson are more relevant phenomenolog- 
ically than its partial decay widths, especially in the intermediate-mass range, where the 
latter cannot be measured directly. In the branching ratios, the universal parts of the 
radiative corrections cancel and, in fact, the only relevant correction is the one specific to 
the H — ► bb decay. 



As a by-product of our analysis, we have confirmed a previous calculation [|46| of the 
0(X t ) correction to the H — > gg decay width. 
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